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A REAPPRAISAL OF THE SOLAR PHOTOSPHERIC C/O RATIO 
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ABSTRACT 

Accurate determination of photospheric solar abundances requires detailed modeling of the solar gran- 
ulation and accounting for departures from local thermodynamical equilibrium (LTE) . We argue that the 
forbidden C I line at 8727 A is largely immune to departures from LTE, and can be realistically modeled 
using LTE radiative transfer in a time-dependent three-dimensional simulation of solar surface convec- 
tion. We analyze the [C I] line in the solar flux spectrum to derive the abundance log e(C) = 8.39 ± 0.04 
dex. Combining this result with our parallel analysis of the [O I] 6300 A line, we find C/0= 0.50±0.07, in 
agreement with the ratios measured in the solar corona from gamma-ray spectroscopy and solar energetic 
particles. 

Subject headings: convection — line: formation — Sun: abundances — Sun: photosphere 



In this paper, we reconsider the carbon abundance pro- 
vided by the forbidden carbon line at 8727 A. Our reanal- 
ysis follows our successful study of the forbidden oxygen 
line at 6300 A, which gave the abundance loge(O) = 8.69 
± 0.05 with a 3D hydrodynamical model photosphere sim- 
ulating the solar granulation (Allende Prieto, Lambert, & 
Asplund 2001a). The principal advantages of a forbid- 
den line over permitted lines are twofold: (i) the line's 
g/-value is most probably more accurately known than 
the ^/-values of the permitted lines, and (ii) the line is 
formed very close to LTE, which is not true for the per- 
mitted lines. We argue that the possible downsides with 
the [C I] line, namely its relative weakness and possible 
blends, can be adequately addressed. 

2. MODEL ATMOSPHERES AND LINE SYNTHESIS 

Our line profile calculations follow very closely those of 
Allende Prieto et al. (2001a) for the [OI] line at 6300 
A. We used a three-dimensional time-dependent hydrody- 
namical simulation of the solar surface as a model atmo- 
sphere (Asplund et al. 2000a and references therein). The 
calculated flux profiles for the forbidden [C I] line at 8727.1 
A and a Si I line at 8728 A result from the average of a 
time sequence of flux profiles computed from 100 snap- 
shots, which are equally spaced over 50 minutes of solar 
time. For each snapshot, the integration over the solar disk 
makes use of intensity profiles for 4x4 angles, accounting 
for the (solid body vsini = 1.9 km s _1 ) rotational broad- 
ening (Dravins & Nordlund 1990). The Uppsala opacity 
package (Gustafsson et al. 1975 with subsequent updates) 
was the source of the continuum opacities, partition func- 
tions, ionization potentials, and other basic data for the 
line synthesis as well as for the simulation. Collisional 

1 We note that Holweger defines the granulation corrections as the difference between the 2D result and the case for the ID spatial average 
of the 2D model atmosphere. Holweger's approach therefore only accounts for the temperature inhomogeneities but not the different overall 
temperature structures in hydrodynamical and hydrostatic models. Finally, we note that there are systematic differences between 2D and 3D 
(Asplund et al. 2000b) 
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1. INTRODUCTION 

Pursuit of quantitative stellar spectroscopy not infre- 
quently introduces new puzzles as old puzzles are probed. 
Recent determinations of the solar photospheric abun- 
dances of carbon and oxygen exemplify this claim. In a 
thorough review of light element photospheric abundances, 
Holweger (2001) recommends the abundances loge(C) = 
8.592 ±0.108, and loge(O) = 8.736 ± 0.078 on the usual 
scale with loge(H) = 12.0. These results based on ID 
NLTE analyses of permitted lines of C I and O I include 
corrections for the effects of solar granulation based on 2D 
modeling 1 . Prior to Holweger's reassessment, the accepted 
carbon abundance was lower (loge(C) = 8.52 ± 0.06) and 
that of oxygen higher (loge(O) = 8.83 ± 0.06; Grevesse 
& Sauval 1998). With the reassessment, the carbon-to- 
oxygen ratio is increased. More intriguingly, the errors 
assigned to the new solar abundances would formally al- 
low a carbon-to-oxygen ratio greater than one. Such a 
ratio is exceedingly puzzling in that it implies remarkable 
abundance differences between the photosphere and the 
corona, the outer planets, and unevolved stars and inter- 
stellar gas in the solar neighborhood. Lowering the oxy- 
gen abundance and keeping the high value for the carbon 
abundance leads also to a conflict with the strength of ob- 
served water and methane bands in the spectra of very 
cool stars (Tsuji 2002). Recognizing the controlling influ- 
ence of the CO molecule on the partial pressures of carbon 
and oxygen, the simple observation of very weak C2 bands 
and strong TiO bands in sunspot spectra shows that the 
carbon-to-oxygen ratio must be less than unity. Yet, a 
precise determination of the photospheric ratio is of great 
interest. 
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Fig. 1. — Comparison between the observed (filled circles) and synthetic line profiles for the [CI] at 8727.1 A. Four parameters have been 
adjusted to minimize the \ 2 '- the wavelength of the [CI] line, a multiplicative correction to the continuum level, the strength of the strong Si 
line at 8728.0 A (log[g/e(Si)]), and the carbon abundance, which is found to be loge(C) = 8.39 ± 0.04 dcx. 



broadening was evaluated using the Unsold's approxima- 
tion (Unsold 1955) for the [CI] line and the neighboring 
Si I line. 

The [CI] 8727.1 A line was identified in the solar spec- 
trum by Lambert & Swings (1967a). They remarked 
upon a potential blending line of Fel, but deemed its 
contribution to be negligible. Kurucz (1993) calculated 
log gf = —3.93 for this iron line, which leads to an equiv- 
alent width of ~ 0.5 mA. Kurucz also calculated the tran- 
sition probabilities for the other 5 lines in the multiplet, 
and four of them appear in the solar spectrum. Using 
the measured equivalent widths (12 < W\ < 40 mA) to 
scale the log g/'s of these four lines, we estimate an equiv- 
alent width of ~ 0.18 ± 0.13 mA 2 . Nave et al. (1994) in a 
thorough reinvestigation of the laboratory Fel spectrum 
did not list the Fe I line with its predicted wavelength of 
8727.130 A, but four stronger lines of the same multiplet 
were measured. Three of the four measured lines appear 
relatively unblended in the solar spectrum. Normalizing 
the line strengths measured by Nave et al. to the relative 
LS-coupling strengths, we can scale the observed equiva- 
lent widths and obtain a third estimate of the equivalent 
width of the 8727.130 A line: 0.11 ± 0.02 mA. Given that 
the solar feature has an equivalent width of about 6 mA, 
the Fe I contribution is very small and can be neglected. 
Gustafsson et al. (1999) eliminate the CN Red system as 
another contributor of a blend. 

The predicted [CI] rest wavelength is 8727.126 A with 
an uncertainty of about 0.015 A (Moore 1993). After cor- 
rection for the gravitational redshift, the center of the fea- 
ture in the solar spectrum is observed at 8727.133 ± 0.009 
A (as determined with the method described in Allende 
Pricto et al. 2002). For weak lines, we expect a signifi- 
cant convective blueshift (Allende Prieto & Garcia Lopez 
1998). Therefore, the predicted wavelength is likely to be 
in error by as much as its full quoted uncertainty. Different 
calculations report results for the Einstein A-value which 



have converged to A = 0.640 s 1 or loggf = —8.136 - see 
Galavis et al. (1997), and Hibbert et al. (1993). A strong 
Si I line affects the continuum in the vicinity of the [C I] 
line. The NIST database quotes a measured wavelength of 
8728.011 A for this transition. We adopted log gf = -0.37 
in order to match the observed line reasonably well. 

3. COMPARISON WITH THE OBSERVED SPECTRUM 

We compare our calculated line profiles with the Fourier 
transform solar flux spectrum of Kurucz et al. (1984). The 
[C I] lies in the scan # 13, which includes the telluric O2 
line at 6883.8335 A, used to normalize the frequency scale. 
This scan has a spectral resolving power in excess of half 
a million, and a maximum signal-to- noise of 2600. 

The Si I line dominates the shape of the local contin- 
uum in the region of the [CI] line. The abundance of 
silicon has been recently studied using the same 3D model 
atmosphere employed here (Asplund 2000). Recent ID 
NLTE calculations suggest that photosphcric Si I lines are 
not much affected by departures from LTE (Wedemeyer 
2001). Since the line's gf- value is uncertain, we consid- 
ered the factor log[g/e(Si)] and a multiplicative correction 
to the local continuum level as free parameters with which 
to adjust the continuum in the vicinity of the [C I] line. A 
variation in the central wavelength of the Si I line has a 
similar effect as a change in the strength of the line. 

The velocity shift between Kitt Peak National Obser- 
vatory and the Sun at the time of the observations has 
been corrected typically within several m s _1 , but other 
factors (setting of frequency scale in FTS, conversion to 
air wavelengths, correction of individual scans with partial 
wavelength coverage, change of velocity shift over long in- 
tegrations, etc.) can degrade the absolute accuracy of the 
wavelength scale up to 0.1 km s _1 (Kurucz et al. 1984). 
The gravitational redshift for photosphcric photons inter- 
cepted at Earth is V g ~ GM /(R c) ~ 0.6336 km s _1 . 
Seasonal variations produced by the ellipticity of Earth's 



2 The quoted uncertainties in the paper represent la values. 
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Fig. 2. — Derived carbon abundance and minimum value of the reduced x 2 as a function of the displacement of the central wavelength for 
the [C I] transition. The best fit is achieved for A[c i] = 8727.139 ± 0.004 A, which is adopted as the zero wavelength shift in the Figure. The 
solid lines are least-squares second-order polynomials fit to the data. 



orbit are about 0.1 m s _1 , Earth's surface gravity induces 
a correction of about 0.2 m s _1 , and uncertainties in the 
solar radius and GMq can hamper V g by about 0.3 m s _1 . 
We are using a hydrodynamical model atmosphere which 
has been shown to predict the convective shifts of pho- 
tospheric lines of weak and moderate strength within 0.1 
km s _1 (Asplund et al. 2000a). Therefore, we should be 
able to predict the measured wavelength of the [C I] with 
an uncertainty of about 0.15 km s _1 (~ 0.004 A), an es- 
timate significantly smaller than for Moore's wavelength. 
We consequently adjust the wavelength of the line to best 
match the observations. 

We proceed to compare our calculated profiles to the 
solar observations adopting different values for the wave- 
length of the forbidden line A[ci]- For each adopted A[ci], 
we minimize the \ 2 between the observed and calculated 
spectra by changing the abundance of carbon loge(C), 
a multiplicative correction factor to the local continuum, 
and log[g/e(Si)]. There are obvious blends on both sides of 
the [C I] line. The line at 8726.8 A is a CN line (Gustafsson 
et al. 1999), but other lines are unidentified. We believe 
there are no blends within the [C I] profile, but, if we are in 
error, the determined carbon abundance is an upper limit. 

First, we fitted 1 A around the [CI] feature, between 
8726.6 and 8727.6 A in the solar atlas, using the Nelder- 
Mead simplex method (Nelder & Mead 1965; Press et al. 
1986). We obtain the best fit, illustrated in Fig. 1, for 
loge(C) = 8.39 dex. We emphasize that the fit is achieved 
without invoking micro- or macroturbulence; the line pro- 
file is predicted from the convective flows without addi- 
tional parameters. By selecting a narrow interval exclud- 



ing the core of the Si I blending line, we are in effect ad- 
justing the line's damping constant to get the optimum 
fit to the wing around the [C I] line; the line depth in a 
damping wing scales as the product of the abundance and 
the damping constant. Next, we repeated the experiment 
for a more restricted wavelength interval, between 8726.9 
and 8727.25 A. The carbon abundance changed by less 
than 0.01 dex, the multiplicative factor to adjust the con- 
tinuum level changed by only 0.0003, and the best values 
for the wavelength of the [C I] by less than 0.001 A. Our 
best fit yielded a reduced \ 2 (37 frequencies and 4 degrees 
of freedom) of 0.87, and, therefore, the chance probability 
for this is P = 0.3. These figures are based on the signal- 
to- noise ratio of <~ 2100 derived by examination of an a; 
parently clean segment of the spectrum around 8731.7 
Fig. 2 shows the variation of the reduced \ 2 an d the car- 
bon abundance as a function of the central wavelength for 
the the [C I] transition, when the more restricted wave- 
length interval is considered. 

The rest wavelength of the [C I] line that minimizes the 
X 2 is 8727.139 ± 0.004 A, after correction for the gravita- 
tional rcdshift. The synthetic spectra predict a blucshift 
of about 0.25 km s" 1 for the line, i.e., a laboratory wave- 
length of 8727.146 A. This differs from Moore's wavelength 
of 8727.126 A obtained not from laboratory measurements 
of the line but from ultraviolet lines connecting the ground 
configuration levels to excited levels. In the introductory 
remarks to her table of C I energy levels, she draws atten- 
tion to Kaufman & Ward's (1966) revisions of ultraviolet 
wavelengths which imply revised energies for the ground 
configuration. These revisions predict a wavelength of 
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8727.141 A in good agreement with our solar-based wave- 
length, even though the uncertainty remains at about 0.01 

A. 

The internal uncertainty in the fit for the carbon abun- 
dance can be estimated as 0.01 dex by confronting the 
expected accuracy in A[crj (0.004 A, which includes the 
observational share of the error), with Fig. 2. Judging 
from the scatter in the most recently calculated A— values, 
the derived carbon abundance is hardly affected by this 
factor (er[loge(C)] ~ 0.005 dex). Other systematic errors 
may arise from the uncertainties in the adopted continuum 
opacities, equation of state, etc., probably about 0.02 dex. 
Finally, if the contribution of the Fe I line discussed in §2 
to the observed absorption is close to the lower limit of 
our estimates, the effect on the derived carbon abundance 
is truly negligible, but if it is close to the upper limit, that 
could decrease loge(C) by ~ 0.03 dex. Conservatively, we 
derive loge(C) ~ 8.39±0.04 dex. As shown by Stiircnburg 
& Holweger (1990), the forbidden line at 8727.1 A is im- 
mune to departures from LTE, and therefore our LTE cal- 
culation of the line formation in a 3D LTE model is likely 
to provide a reliable abundance. 

It is interesting to compare our 3D result to a stan- 
dard analysis using one-dimensional model atmospheres. 
Use of the solar semi-empirical model atmosphere derived 
by Allendc Prieto et al. (2001b), and the micro- and 
macro-turbulence derived with the model (£ = 1.1 km 
s _1 , macro=1.54 km s _1 ), leads to loge(C) = 8.47 dex, 
the Holweger-Miiller model gives loge(C) = 8.48 dex, and 
a flux-constant MARCS model atmosphere loge(C) = 8.41 
dex. 

4. THE SOLAR CARBON ABUNDANCE AND C/O RATIO 

Our determination of the solar carbon abundance 
(loge(C) = 8.39) from the forbidden line at 8727.1 A is 
lower than the most recent estimates based on this line. 
Our lower abundance is mainly caused for two reasons: a 
revised /—value (~ —0.07 dex), and the use of a 3D model 
atmosphere (~ —0.08 dex). 

Other indicators of the solar photosphcric carbon abun- 
dance are available. CI, CH, and C2 lines have received 
significant attention in the literature and suggest, when 
analyzed with empirical or theoretical ID model atmo- 
spheres, higher abundances than our determination from 
the [C I] line (see, e.g. Grevesse et al. 1991). Pending 
a full reanalysis with the 3D model, we note that the C 
abundance from the molecular lines will be reduced in the 
change from a ID to 3D model due to the high temper- 
ature sensitivity of molecule formation and the presence 
of cooler gas in the 3D model (Asplund & Garcia Perez 
2001). To assess the effect on the CI lines, we selected 
the weakest and least blended 13 lines in overlapping lists 
of solar lines used previously by Biemont et al. (1993) 
and Stiircnburg & Holweger (1990, also Holweger 2001). 
Analysis of the published equivalent widths for the cen- 
ter of the solar disk spectrum with the ^/-values in the 
NIST database and the 3D model gives the LTE abun- 
dance loge(C) = 8.44 ±0.06. NLTE 3D calculations are 
not available, but if the ID corrections for NLTE effects 
from Stiircnburg & Holweger are adopted, the abundance 
is lowered to loge(C) — 8.42 ± 0.06, a value in harmony 
with our result from the [C I] line. Our preliminary conclu- 



sion is that the lower carbon abundance is not contradicted 
by other indicators. 

Combining the updated solar abundances of carbon and 
oxygen from the forbidden lines, we obtain a number- 
density ratio C/0= 0.50 ±0.07. Independent estimates of 
the C/O ratio can be obtained for the solar wind and the 
corona. The C/O ratio has been measured from gamma- 
rays produced by the interaction of particles produced in 
solar flares with the surrounding medium. Solar Maximum 
Mission (SMM) data for a flare observed in 1981 yielded 
C/O = 0.42 ± 0.09 (Murphy et al. 1990; Fludra et al. 
1998). Ramaty et al. (1996) derived 0.35 < C/O < 0.44 
for 19 events observed by SMM. More recently Murphy 
et al. (1997) obtained 0.54 ± 0.04 from CGRO data of a 
flare observed in 1991. Solar energetic particles detected 
near Earth can themselves be studied to determine relative 
abundances. Reames (1995) reviews coronal abundances 
derived from solar energetic particles. Compared to pho- 
tosphcric abundances, elements with a low first-ionization 
potential (FIP) (<; 9 eV) are enhanced in the corona, but 
those with a higher FIP arc not. As carbon and oxygen 
have a similar, high, FIP, this effect is not expected to dis- 
turb their ratio greatly from its photosphcric value. For 
'gradual' events associated with coronal mass ejections, 
the C/O ratio shows very little variation from one event 
to another: Reames gives C/O = 0.465 ± 0.009. For events 
associated with flares, C/O = 0.434 ± 0.030. These results 
are in good agreement with our 'forbidden' photospheric 
value. 

The revised abundances of carbon and oxygen in the 
Sun compare very well with the abundances (gas plus 
dust) derived from recombination lines in the Orion neb- 
ula (Esteban et al. 1998): loge(C) = 8.49 ± 0.12, and 
loge(O) = 8.72 ± 0.07 (C/O = 0.59). They are also in 
good agreement with the average values found in B-type 
stars in the field: loge(C) = 8.31, and loge(O) = 8.58 
(C/O = 0.54), as well as in clusters: loge(C) = 8.35, and 
loge(O) = 8.69 (C/O = 0.46) (Kilian 1992, 1994; Adel- 
man Robinson & Wahlgrcn 1993; Gies & Lambert 1992; 
compiled by Snow & Witt 1996). 

Revision of the solar C and O abundances using their 
forbidden lines brings these abundances into line with 
recent results for B stars and the Orion nebula. One 
may speculate that this close correspondence in compo- 
sitions may extend to other elements, for example, nitro- 
gen, through revisions to either the solar, stellar or nebular 
abundances. Unfortunately, the [N I] lines are too weak for 
positive identification in the solar spectrum (Lambert & 
Swings 1967b). Nitrogen is represented by weak NI lines 
for which Holweger (2001) gives loge(N) = 7.93 ± 0.11 
after including NLTE effects and a correction for granula- 
tion. Combining the N/O ratios measured for solar ener- 
getic particles with the photosphcric abundance of oxygen 
gives a nitrogen abundance which is about 0.1 dex lower 
than Holweger's, and in better agreement with nearby B 
stars and the gas in Orion. 
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